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O I We present the Spitzer Infrared Spectrograph (IRS) spectrum of SR20, a 5-10 

I AU binary T Tauri system in the p Ophiuchi star forming region. The spectrum 

has features consistent with the presence of a disk; however, the continuum slope 



is steeper than the A slope of an infinite geometrically thin, optically thick 



Oh' disk, indicating that the disk is outwardly truncated. Comparison with photom- 

O ' etry from the literature shows a large increase in the mid-infrared flux from 1993 

I to 1996. We model the spectral energy distribution and IRS spectrum with a wall 

^ ' + optically thick irradiated disk, yielding an outer radius of 0.39],]o;q3 AU, much 

^ smaller than predicted by models of binary orbits. Using a two temperature 

>- ■ minimization model to fit the dust composition of the IRS spectrum, we find 

^ ■ the disk has experienced significant grain growth: its spectrum is well-fit using 

' opacities of grains larger than 1 /im. We conclude that the system experienced 

a significant gravitational perturbation in the 1990s. 
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Introduction 



Classical T Tauri stars are typically surrounded by optically thick accretion disks, indi- 
cated by strong and broad Ha emission lines and by a characteristic infrared excess. Over 
the last twenty years, observations have revealed that many of these young; stars are com- 



pone nts in binary, and sometimes higher order, systems (e.g. iSimon et al.lll992l : iGhez et al. 
19931 ) . A binary system may pr oduce up to three distinct disks: a circumprimary, cir cumsec- 
ondary, and circumbinary disk (lArtymowicz fc Lubowll 19941 : 1 Jensen fc Mathieulll997l . (AL94; 
JM97)). Gravitational interactions between binaries affect the structure of their disks; in 
particular, the maximum outer radius for a circumstellar disk in a binary system is approx- 
imately 0.18 to 0.4 times the semi-major axis (AL94), depending on the mass ratio and 
orbital eccentricity. 

SR20 is one such close binary in the p Ophiuchi star forming region. I t s com ponents 
have an angular separation ranging fro m 0.038 to 0.071^^f[Ghez et al. 1993 . 1995 ) which, 
at a distance to Ophiuchus of 140 pc (Ide Zeeuw et al.l Il999l ) corresponds to a projected 
linear separ ation of 5.3-9.9 A U. The secondary is 2.2 magnitudes fainter than the primary 
at 2.2 /xm (IGhez et al.l Il993l ) and therefore is less massive than the primary, diskless, or 
both. Ha emission ha s been detected from SR20 w i th equivalent width s ranging from 15 



to 21 A( Rydgren 1980 : Bouvier fc Appenzellerl 



1992 



Wilking et al.ll2005l . (W05)). Further 



observations by iTakami. Bailey fc Chrysostomoul (120031 ) did not find a spectro-astrometric 
signature in the Ha line, indicating that there is little to no accretion in the secondary relative 
to the primary. Based on the lack of evidence for a circumsecondary disk, we assume that 
the infrared excess of SR20 is dominated by a circumprimary accretion disk. 



Here we present Spitzer Space Telescope Infrared Spectrograph (IRS) (IHouck et al. 



20041 ) observations of the SR20 system. We model the structure of the system and the 
dust it contains and analyze the nature of the SR20 system. While instances of sculpture of 
proto-planetary accretion disks b y planetary or stell ar companions within the disk are be- 
coming widely known (e.g. JM97. iFurlan et al.ll2007l ). SR20 is an apparently rarer example 
of sculpture from without. 



2. Observations and Data Reduction 

We observed SR20 on 2006 April 16 (Spitzer AOR 12698368), with the short wavelength 
(SL; 5.2-14 /im) and long wavelength (LL; 14.0-36.1 iim) low spectral resolution (A/AA=60- 
120) IRS modules. At each order of each module, the object was observed twice, once in each 
of the nominal nod positions, 1/3 of the way from the ends of the slit, with exposure times 



3 



for each observation of 6 seconds in SL and 14 seconds in LL. The spectrum was extracted 
from the Spitz er Science Center (S SC) S14.0 pipehne basic cahbrat ed data using t he SM ART 
software tool (IHigdon et aLll2004l ) following the same procedure as lFurlan et al.l (120061 ) with 
the exception of background subtraction, which was done using the opposite nod position. 

Assuming nod-to-nod flux differences were caused by slight mispointing, for each module 
we scaled the nod with the lower flux to match the flux level of the nod with the higher 
flux. Scale factors were typically within a few % of unity. The resulting nods were averaged 
to obtain the final spectrum, and the spectral uncertainties are estimated to be half the 
difference between the two independent spectra from each nod position. We estimate the 
spectrophotometric accuracy of the final result to be around 5%. 



3. Analysis 



3.1. Spectral Energy Distribution 

The 5-36 /im IRS spectrum of SR20 (Figure [T]) has a clear infrared excess and weak 
silicate dust emission at 9.8 and 18.0 /im, consistent with the presence of a circumstellar 
disk. This excess cannot be attributed to an optically thin shell, as in mass-losing giant 
stars, since the 5-8 fim H2O absorption features indicate that the underlying continuum is 
optically thick. 

Here we adopt a spectral type of G7 and Ay = 5.7 mag, determined from optical data by 
W05. To correct the photometry and IRS spectrum for continuum slope fittin g, we applied 
a "standard" extinction correction model with Ay/rg.y = 12.5 and Ry = 3.1 (jDraind l2003l . 



D03). For the dust composition fitting, we corrected more carefully over the silicate features 
by scaling our previous correction model to Ay/Tg^ = 25, which may be more appropriate 
for silicate features in dark clouds (IChiar et al.ll2007l ). For consistency, we compared results 
for both dust composition and disk structure and they are independent of the correction 
used. 

After exinction correction, we constructed two SEDs using the Ay/rg y = 12.5 corrected 
spectrum. The first contains optical data, the near infrared data closest in time to our 
spectru m, the IRS spectrum and a 5600 K, solar metallicity NextGen stellar photosphere 
model (IHauschildt et al.lll999l ) (Figure [2]). We chose to normalize the photosphere model to 
the 2MASS J-band flux density rather than that at I-band due to the greater variability at 
I (1.0 magnitude vs. 0.4 magnitudes at J). It is unlikely that SR20 suffers from significant 
veiling at J, since the line strengths in the iLuhman &: Riekd (119991 ) K-band spectrum are 
comparable with a late G photosphere. The second SED shows the near infrared data, sorted 
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into four epoch "bins": 1974, 1993, 1996-1997, and 1999-2004 (Figure [3]). Significantly, in 
1974 and 1993 the excess from 2 to 10 fim was small but increased substantially by 1996. 
Changes in the 2-10 fim excess since then have been small. 

It can be seen from Figure 2 that SR20 exhibits excess emission above the photosphere 
over A = 1.25-36 and that the slope of the continuum emission is much st eeper than 
A^^, the slope of an infinite geometrically-thin, optically-thick irradiated disk ( Hartmann 
19981 ). There is a correlation between the continuum slope and dust settling in the disk; for 
the IRS Taurus sample, the bluest and presu mably most settled T Tauri disks had continuum 
slopes approaching, but not equal to, A~^/^ (IFurlan et al.ll2006l ). as demonstrated by the IRS 
spectrum of IS Tau, also in Figure [2l We interpret a slope steeper than A~^/^ to indicate 
that the circumstellar disk is outwardly truncated. 



3.2. Wall + Disk Model 

To investigate whether outward truncation is plausible and determine the degree to 
which it occurs, we modeled the disk over a range of radii. We began by modeli ng the inner 



edge of the disk with a vertical "wall" loca ted at the dust destru ction radius (INatta et al. 



200ll ). calculated following the methods of iD'Alessio et al.l (120051 ) . which include the wall 



atmosphere. The wall emission depends on its radius, Rin, height Z, and the inclination to 
the line of sight, i. In turn, Rin depends on the dust composition and size distribution in 
the wall and on the sum of the stellar, L*, and accretion. Lace, luminosities. We use a size 
distribution oc a~^-^ between a^i^ = 0.005/im and amax- For L*, we adjusted th e value from 



W05 to 140 pc, yielding 6.76 Lq. For a list of stellar parameters, see Table [TlNatta et al. 



( l2006l ) estimated an upper limit to Lace of 0.03 Lq, using the luminosity of the Pa/3 line. 
They adopted an extinction of Aj = 1.3, which corresponds to Ay = 4.6, lower than the 
one used here. Scaling up the J band flux according to the difference of extinction values, 
we get an upper limit to the accretion luminosity of 0.05 Lq, which is much lower than 
L*. Therefore, we neglect the accretion luminosity in the calculation of the dust destruction 
radius and use as input parameters the maximum grain size, amax, the inclination, and the 
height of the wall. Varying amax changed Rin and the peak-to-continuum ratio of the silicate 
features, which formed in the wall atmosphere, with large Omax producing a wall radius 
closer to the star and a lower peak-to-continuum ratio. We were able to fit most of the SED 
with just the emission from a wall with spherical solid amorphous olivine grains that have 
grown to 1.0 fim. With these grains, the dust destruction radius is locat ed at 0.34 AU. In 



the v ertical wall approximation, the maximum emission occurs for ~ 50° (jPuUemond et al. 



20011 ). At this inclination, best fit is for a wall height Z = 2.11 R,,, or 2.6 scale heights. 
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which is lower than the expected range of 4-5 scale heights (iDullemond et al.ll200ll ). Lower 
inclinations would require higher walls to fit the near-infrared excess. The wall emission 
alone is insufficient to explain the IRS spec trum past 10.0 Mm , so we added a very small fiat, 
optically thick irradiated disk component (IHartmannl Il998l ) . The temperature distribution 
of the disk component is given by 



T(r) 



lo 



sin^^^ - ^ f 1 - ^ 



where r is the cylindrical radius, and Iq 



. From 1.25-971 /xm, we numerically integrated 



the corresponding fiux density over the range of radii from i?in, the wall radius, to -Rout, the 
outer radius. By adding a disk component with R^^ = 0.34 AU, i?out = 0.39^q;o3 AU, and i 
= 50° to the wall model at the same inclination with 1.0 /xm solid olivine grains, the SED 
of SR20 is well fit (Figure H]). The only discrepancy is on the longer wavelength side of the 
silicate features, which could be fit better with slightly larger grains. As noted, we could get 
similar fits in the near-IR with lower inclinations and taller walls. However, the emission of 
the optically thick disk scales as cost, increasing as the inclination decreases. As a result, to 
obtain the same fit to the mid-infrared, the disk radius would have to be smaller. A similar 
conclusion could be expected if the w all were round and its maximum emission occurred at 
low inclinations (jlsella fc Nattall2005l ). Therefore, 0.39 AU is an upper limit for the radius 
of the disk. 

To understand the relationship between the binary companion and circumprimary disk, 
we need information about the binary orbit. We were unable to determine a unique solution 
(there are only three separation measurements in the literature), but we placed rough con- 
straints on the semi- major axis and eccentricity. In the limit of the eccentricity approaching 
unity, the semi-major axis is half the maximum separation. Therefore, regardless of the 
eccentricity, the semi-major axis of SR20 must be greater than 5 AU. For our preliminary 
orbit estimate, values between 8 and 12 AU fit best along with lower eccentricities of 0.2-0.3. 



3.3. Dust Composition Model 



Although our wall + irradiated disk model produced a good fit to the IRS spectrum, 
it uses only one silicate component at a time with a fixed set of Omax- To perform a more 
detailed analysis, we calculate a two temperature model of the silicate dust composition by 
minimization with respect to weights proportio nal to dust mass, bl a ck body solid angle . 



and temperature for black bodies and dust gr ains (Sargent et a 



Chen et al.l l2006l ) with optical constants from iDayl (1l979l ) and 



20061: iKastner et al.ll2006 



Koike et al 



(|2003h . 



Grains 
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larger than 8 /im in radius, which do not contribute to the sihcate features, could still 
contribute to the black body components. The best fit {xu = 4.0) occurs at 410 and 1221 
K with the components listed in Table [H which are shown in Figure [TJ Interestingly, the 
only significant components besides the black bodies are large (5 /im) porous (60% vacuum), 
amorphous olivine (Mg2Si04) grains, which is consistent with the 10 fim peak to continuum 
ratio of 1.35:1. Submicron amporphous olivine and pyroxene, crystalline silica, and large 
amorphous pyroxene, all of which are typical of T Tauri disks, were included in the modeling 
but found not to be present. Forsterite was present in the best-fit model, but at less than la 
levels, so we do not consider it to be significant. For d = 140 pc, our fit yields 9.55 x 10^^ 
^lunar of large amorphous olivine grains, combining the contributions at both temperatures, 
with an upper limit of 1.36 x 10~^ Miunar on sub-micron grains. 



Discussion and Conclusions 



The SED, system geometry, and dust composition shed light on the nature of the SR20 
system, while simultaneously raising further questions. 

The IRS spectrum is well-fit by our composite wall -|- circumprimary disk truncated 
beyond 0.39 AU. Outward truncation could be caused by the orbit of a companion or lack 
of replenishment from a circumbinary disk. Using the tidal truncation models of AL94, the 
range of likely eccentricities (0.2-0.3) and semi-major axes (8-12 AU) of the A-B orbit are 
not extreme enough to cause truncation. Even with the lower limit on the semi-major axis, 5 
AU, and the highest order resonances and eccentricities tested in AL94, the minimum stable 
outer truncation radius caused by the known binary is 0.9 AU, considerably larger than the 
0.39 AU we derive (AL94). However, the AL94 models may not be appropriate for a disk 
with such a small ratio between the outer and inner radii, so we cannot absolutely rule out 
truncation by the known binary, at least on this basis. Alternatively, it may be that SR20 
is a hierarchical triple. A companion object with an eccentricity < 0.8 could orbit in the 
region between 0.98 AU and 2.17 AU from the primary, causing a disk truncation at 0.39 
AU. Such a close, low mass companion would be difficult to detect. 

Based on the 0.39 AU truncation radius, we speculate that replenishment of the cir- 
cumprimary disk by a circumbinary disk occurred in the past. A disk betwe en 0.1 and 10 
AU would accrete on a timescale of 1,000-100,000 years (jQuillen et al.ll2004j ) and the me- 
dian age of the Ophiuchus cluster is 2.1 Myr (W05). SR20 has been observed at 800, 850, 
and 1300 fim with only 3 a up per limit results, implying less than 5.4 x 10^ Miunar (540 
Miunar of dust) in disk material ([Andrews fc Williams! 120071 ). Low cir cumbinary disk ma sses 
are typical of binaries with semi-major axes between 1 and 50 AU (I Jensen et al.lll994l ). so 
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submillimeter non-detection does not indicate a depleted circumbinary disk; however, it is 
unlikely that lack of replenishment is responsible for the 0.39 AU truncation radius. 

A nearby companion truncating the circumprimary disk is perhaps also the best ex- 
planation for the photometric variability. Before 1993, the infrared excess in SR20 started 
in the N band, indicating that there were few small dust grains close to the central star. 
However, the Ha emission indicates that a gaseous accretion disk, perhaps mixed with much 
larger planetesimals, was present. Between 1993 and 1996, the infrared excess increased 
sharply, indicating an increase in the amount of dust in the disk, which could be explained 
by a gravitational interaction between the known binary and a third, unseen, body, causing 
the unseen companion to orbit closer to the primary, perturbing the circumprimary disk. 
Coupling between the gas and dust could confine the newly created dust to the accretion 
disk. The appeareance of the observed excess over 3 years, preceded by 20 years without an 
excess and followed by 10 years of little variation in the observed excess, seems consistent 
with truncation by a nearby source with a shorter period than is plausible for the known 
companion. 

The silicate dust composition of the system is not typical of a Class H object; our fit 
indicates that the optically thin portion of the disk is comprised primarily of large grains 
(> 1 fim). Although having a large fraction of large grains is consistent with debris disks o r 



advanced dust processing in the inner regions of a protoplanetary disk (IChen et al.ll2006l ). 
this is the first time we have seen a Class II disk with negligible submicron grains. Taken 
with the aforementioned explanations of the 1993-1996 variability, the large grains could 
support a collision scenario: gravitional interactions between the known secondary and the 
unseen companion cause collisions between planetesimals, producing dust grains that are 
detected as the current disk. If there were a population of planetesimals - invisible in our 
spectra - comprising an extension of the disk beyond the outer truncation radius of the dust, 
truncation by the known companion would be more probable. 

Near infrared interferometry and orbital monitoring would be necessary to determine 
the nature of this complex system more precisely. 



This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology, under NASA 
contract 1407 and made use of the ADS, SIMBAD, and Vizier utilities. Support was provided 
by NASA through contract number 1257184 issued by JPL/Caltech, JPL contract 960803 
to Cornell University, and Cornell subcontracts 31419-5714 to the University of Rochester. 
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Wovelength (/xm) 



Fig. 1. — Dust composition model overplotted on the extinction corrected IRS spectrum. 
Light grey and dark grey represent 410 K and 1221 K components respectively. Compo- 
nents are black bodies (dash/dot), large amorphous olivine (solid), and forsterite (dotted). 
Forsterite is present at less than la strength and not considered significant. 
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Fig. 2.— "Current" SED of SR20, with 1983-84 UBVRI jBouvier fc Appenzelleil [l992l ). 
1999 JHK (2MASS), 2004 IRAC (C2D), and MIPS 24 /xm photometry with our 2006 IRS 
spectrum. Data were extinction corrected using Av=5.7 and the DOS law. UBVRI data are 
plotted with bars to indicate the extent of variability. Overplotted are the IRS spectrum of 
IS Tau (x 15) and a line representing A"^/^. 
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Fig. 3. — SR20 photometry from four epochs: 1974: JHKL (iRydgren. Strom fc Stromlll976l ): 
1993: J HKL'M'N fJM97j: 199 6-7: JK (DENIS), 3.6 /xm (ISO-PHT, sky+source), 6.7 and 
14.3 /im JSontemps erZlbnOlh : 1999-2004: JHK (2MASS), IRAC (C2D), 24 /xm (Extracted 
from MIPS AOR 4321280 with an apertur e of radius 14.94 " sky annulus from 29.88 to 42.33" 
and aper ture correction factor of 1.143 ( Su et al. 20061 ) using the SSC's APEX software 
package (IMakovoz fc Marleaull2005l )). Photometry is plotted with error bars, most of which 
are smaller than the plot symbols. 
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Fig. 4.— Best fitting model (grey sofid fine): d = 140 pc, i = 50° i?in = 0.34 AU, and R, 
= 0.39 AU. 
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Table 1. Stellar and model properties 



Component 


Value 


Reference 


Stellar parameti 


srs 




SpT 


G7 


W05 


Te (K) 


5584 


W05 


i. (L©) 


6.76 


W05, adjusted" 


(Mq) 


1.86 


Siess et al. (2000) 


R* (R©) 


2.57 


Sess et al. (^2000) 


Lace {L© 


< 0.05 


Natta et al. ("2006). adlusted" 


Ay (mag) 


5.7 


W05 


Dust Components 


> Icr 




Parameter 


Value'' 


1(7 uncertainty'^ 


(410 K) 






black body 


7.12 X 10-1'' 


3.84 X 10-1* 


large Mg2Si04'' 


3.41 X 10-20 


5.18 X 10-21 


(1221 K) 


black body 


2.56 X 10^1' 


4.70 X 10- 1'-* 


largo Mg2Si04'' 


3.58 X 10-21 


5.97 X 10-22 



Note. — " See section 3.2 for details. ^ Blackbody components are solid angles in units of sterradians, while dust components 
are column density times solid angle in g/cm2. Uncertainties are the increments away from the best fit values enough to 
increase Xf t>y 1.0. Amorphous olivine 



